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ABSTRACT 

We explore the physics of crystallization in the deep interiors of white dwarf stars using the color-magnitude 
diagram and luminosity function constructed from proper motion cleaned Hubble Space Telescope photometry 
of the globular cluster NGC 6397. We demonstrate that the data are consistent with the theory of crystallization 
of the ions in the interior of white dwarf stars and provide the first empirical evidence t hat the phase tra nsition 
is first order: latent heat is released in the process of crystallization as predicted by Ivan H orn (1968)). We 
outline how this data can be used to observationally constrain the value of Y = ^Coulomb /^thermal near the onset 
of crystallization, the central carbon/oxygen abundance, and the importance of phase separation. 
Subject headings: white dwarfs — dense matter — equation of state 



1 . STAR FORMATION HISTORY AND PHYSICS FROM THE WHITE 
DWARF STARS 

White dwarf stars are the inevitable progeny of nearly 
all (~ 97%) stars (e.g jFontaine. Brassard. & Bergero nl200lL 
hereafter FBB). Their distribution can be used to extract two 
things: age of the stellar population and cooling physics of the 
white dwarf (WD) stars. The two are interrelated, but quali- 
tatively different. Extracting the age and history of star for- 
mation has become known as WD cosmochronology. An ex- 
cellent review emphasizing this connection and the attendant 
uncertainties is given by FBB. 

The techniques of WD cosmochronology have been suc- 
cessfully applied to the disk by a number of investigators ( e.g., 
IWinget et alJIl 9871: IWoodfl 9921: lHansen & Lieberfeool and 
references therein) and are being continuously refined. They 
have also been applied to a variety of open clusters and cal- 
ibrate d against main sequence turnoff and relate d methods 
(e.g.. iKalirai et all 120071 iDeGennaro et all 12008, and refer- 
ences therein). The Hubble Space Telescope (HST) photom- 
etry obtained by Richer and Hansen and thei r collaborators 
(lHansen et alJl2002t l2007t iRicheret al.ll2008l) has yielded a 
new harvest of information for WD populations. They have 
used the Advanced Camera on the HST to reach the terminus 
of the WD cooling sequence, giving us a quali tatively dif- 
ferent tool for analyzing the WD population. lHansen et all 
(2007) used Monte Carlo techniques in conjunction with their 
cooling models to determine the age of NGC 6397 from the 
WD stars, attempting to account for uncertainties in the basic 
physical parameters of the WD stars to determine an age for 
the cluster; using goodness of fit criteria, they arrive at an age 
for the cluster, based on WD cooling, of 1 1 .47 ± 0.47 Gyr. 

Finding the signature of the key physical properties of the 
WD stars in the disk luminosity function (hereinafter LF) has 
proven more difficult than getting an age constraint. This is 
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because the disk population contains stars formed at different 
times and from different main sequence progenitors. This is 
greatly simplified in a cluster sample and, most of all, in an 
old globular cluster. In this paper we focus on the HST pho- 
tometry of NGC 6397 and the distribution of WD stars in the 
color magnitude diagram. We report evidence for a "bump" 
in the LF due to the release of the latent heat of crystallization 
and we show how this can be refined to yield more accurate 
measures of these processes. 

2. ANCHORING THE WHITE DWARF EVOLUTIONARY 
SEQUENCES IN THE COLOR-MAGNITUDE PLANE 

We fit main sequence, pre- WD, and WD evolutionary mod- 
els simultaneously. The main sequence and pre-WD models 
we used for this wo rk were computed with the Padova stel- 
lar evolution code (Marigo et al. 2008). We used a variety of 
metallicities to determine the best fit to the main sequence and 
WD models. 

Our WD evolutio nary models have updated constitu- 
tive p hysics (see e.g.. lBischoff-Kim. Montgomery. & W inget 
2008). We place the new generation of WD evolutionary 
models of DA and DB WD stars in the observed Fsuw 
vs. F606W color-magnitude diagram (hereinafter CMD) using 
P. Bergeron 's model atmosphere grids 7 (for a detailed de- 
scription seelBergeron et alTl9 95; Holberg & Bergero nl2006l ; 
Holberg, Bergeron, & Gianninas 2008) along with an analyt- 
ical correction to the Kowalski (2007) results for the effect of 
Ly a far red-wing absorption. This correction is small and 
will be discussed in a forthcoming paper. 

Our best fit of the CMD in the natural ACS color sys- 
tem gives a metalhcity of Z = 0.00012 ±0.00001, £(606W- 
814W) = 0.22 ±0.02, and (m-M) = 12.49 ±0.05 (Figure 1) 
and a main sequence turnoff age of 12^^_Gyr. The a ge is 
consistent wi t h the values found by Richer et al. (2008|) and 
lHansen et al.1 d2007l) . even though they used the Dartmouth 
Evolutionary Se quence (PES). The metallicity is a factor of 
two lower than iRicher et al.1 (120081) but is in agreement with 
the indepe ndent direct spect roscopic determinations based on 
VLT data (Korn et al. 2007). The values of these parameters 
fix the WD cooling tracks in the color-magnitude diagram and 
eliminate the freedo m to slide the track s as has been done in 
other works (e.g., lHansen et alj|2007b . This constrains the 

7 http://www.astro.umontreal.ca/~bergeron/CoolingModels/ 
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Figure 2: NGC 6397 white dwarfs with DA and DB evolutionary tracks with 
Bergeron-Kowalski Colors and E=0.22,{m-MH2.49 



12 3 

F606W-F814W 

FIG. 1. — The best fit to the proper-motion screened HS T data on 
NGC 6397 using the Padova stellar evolution code i Marigo et al. 2008) and 
the Bergeron (2006) atmospheres. This fit gives [Fe/H] = -2.2, (m-M) = 
12.49, and E(F6Q6W - FSUW) = 0.22. This anchors our WD sequences in 
the CMD. 

best-fit total WD mass throughout the CMD. It is clear from 
Figure 2 that the tracks fit the bulk of the sample well. 

The data (see Fi gure 2) are taken fro m the proper motion 
selected sample of iRicher et alj ( 20081). This i s a m ore ho- 
mogeneous sample than that of lHansen et a l. (2007) but is 
smaller because of the reduced area and magnitude limits 
of the proper-motion data. Four features stand out in Fig- 
ure 2: there is a gap in the distribution near F814W = 24.5 
that may be statistically significant, there is a noticeable con- 
centration of stars near F814W = 26.5, there is a terminus 
at approximately F814W = 27.6 and a noticeable turn to the 
bl ue before the terminu s. These last two features were noted 
in lHansen et al.l (120071) . In this paper we focus on what we 
can learn from the concentration, or clump, of stars near 26.5, 
providing a physical explanation. 

3. PHYSICS WITH THE CMD AND LUMINOSITY FUNCTION 

The CMD diagram cons trains the mechani cal and thermal 
properties of the WD stars (Ri cher et al.ll2008l) . Once the evo- 
lutionary tracks have been anchored by the main sequence 
and WD sequence simultaneously as described above we can 
move on to explor ing the physics contained in the CMD. 
lHansen et al.l (12007) point out that the location of the terminus 
provides a simple lower-limit to the age of the cluster from the 
WD cooling times. For our models, this WD cooling limit is 
reached at about 10.5 Gyr for pure carbon core models. This 
age lim it is consistent wi t h the values quoted in lHansen et al.l 
(120071) and iRicher etal] (120081) . The position of the tracks 
is insensitive to processes affecting only the age; to examine 
these we must look to the LF, the number of stars observed as 
a function of magnitude. 

For the parameters described above, it is evident from Fig- 
ure 2 that models with masses in the range of 0. 500-0. 535M Q 
best fit the region near the center of the clump. This increases 
slightly for decreased values of m-M and is also a function 
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FIG. 2. — NGC 6397 WDs with DA and DB evolutionary seque nces using 
Bergeron (2006) atmospheres adjusted with an analytical fit to the Kowalski 
1 2007) changes for Ly-a red-wing opacity. 

of the reddening. We emphasize again that simultaneously 
fitting the main sequence and white dwarf sequence provides 
tight constraints on both the distance modulus and the redden- 
ing 



I n Figure 3 we show the LF of both the Han sen et al.l d2007l) 
and Richer et al.j (I2008I) samples. The peak in both LFs near 
F814W = 26.5 suggests t hat evolution slows through this re- 
gion. In the sample of lHansen et alj (|2007) the LF contin- 



40 - 




24 26 
F814W 



FIG. 3. — top panel: Th e observed W DLF (histogram) and completeness 
relation (dotted line) of thelRicher et alj J2008h sample for NGC 6397. lower 
panel: the same for the Hansen et al. (2007) sample. For both samples, we 
note that the completeness changes fairly slowly over the region of the ob- 
served rapid fall-off of stars and while remaining above 50%. 
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ued to rise well past this peak, so it was not the maximum of 
the distribution; it completely dominates the distri bution after 
the ap plication of the proper motion selection (Rich er et alj 
2008). The completeness estimates of both samples have been 
carefully considered by the respective authors (Figure 3, dot- 
ted lines) and the relatively slow variation of the completeness 
near this peak implies it is not the result of incompleteness. 
We therefore seek a physical explanation of this peak in the 
context of a physical process that occurs near this point in the 
models. 

Two processes occur in the dominant DA models near this 
point: crystallization and the convective coupling (e.g., FBB). 
Crystallization, through the release of latent heat, slows down 
evolution and produces a bump in the LF. Convection, when 
it reaches down to the degeneracy boundary, decreases the in- 
sulation of the nondegenerate envelope and temporarily in- 
creases the total temperature gradient; this serves to slow 
down the evolution, briefly, then causes it to speed up again. 
This produces a broad feature in the WD cooling curve that 
will have a signature in the LF. 

3.1. The r of Crystallization 

Crystallization in the dense Coulomb plasma of WD inte- 
riors w a s theoretica l ly pred icted i ndependently by iKirzhnitsI 
(1960), Abrikosov (I960]), and ISalpeterl |l961) to occur 
when the ratio of the Coulomb energy to the thermal en- 
ergy of the ions (the ratio 'T") is large. For a one- 
component plasma (OCP), there is universal agreement 
among different theoretical approaches that crystallization oc- 
curs when r ~ 175 (e.g.. ISlatterv, Doolen. & Dewm1ll982l, 
String fellow. Dewitt, & Slattervl 119901; iPotekhin & Chabrierl 
2000; Hor owitz. Berry. & Brownl 12007). There is a similar 
consensus, based largely on a density-functional approach, 
that this value of T also holds for a binary carbon and oxygen 
mixture. Such a mixture is likely re l evant to WD interiors. 
Recently, iHorowitz. Berry. & Brownl ((2007) used a massive 
molecular dynamics computation to explore crystallization in 
a dense Coulomb plasma. They found T = 175 for an OCP, 
while for a specific mixture of elements they found T ~ 237. 
As we show, such a difference is potentially measurable from 
the observations of WD stars in globular clusters or older open 
clusters. 

3.2. Luminosity functions with and without crystallization 

lHansen et alJ (|2007) demonstrated in their analysis that the 
entire observed sequence represents a very narrow range of 
WD masses, including magnitudes well below F814W = 26.5. 
It is therefore reasonable for purposes of this initial explo- 
ration to adopt a fiducial mass. On the basis of the model 
tracks in the CMD shown in Figure 2, we choose the model 
that passes nearest the color of the red edge of the clump 
of stars corresponding to the peak in the LF; this model 
has a mass of 0.5M©. For the layer masses, we assume 
M H /M* = 10" 4 and M He /M* = 10" 2 for the DA sequences, 
and Mne/M+ = 10~ 2 for the DB sequence. We adopt a car- 
bon core model including the effects of crystallization for this 
sequence. 

Assuming a constant star formation rate, the theoretical LF 
is proportional to the "cooling function" of an evolutionary 
model sequence. This function is given by the derivative 
(dt I dm), where m is the F814W magnitude of a given model 
and t is its age. Since we will be comparing directly with the 
data, we also multiply the theoretical LF by the completeness 
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FIG. 4.— The observed WDLF of NGC 6397 (Richer et alj|2008l his- 
togram) with LFs from theoretical evolutionary sequences of 0.5Mq DA 
models with pure carbon cores (lines): crystallization with ("Pure_C_Xtal") 
and without ("Pure_C_Xtal_No_LH") the release of latent heat, and exclud- 
ing the physics of crystallization altogether ("Pure_C_No_Xtal"). The nor- 
malization of the theoretical curves is chosen to minimize the RMS residuals 
in the neighborhood of the peak, between the magnitudes of 25.1 and 27.7, 
the faintest value calculated for the no crystallization case. The value of the 
average residual for each curve is listed in the legend, e.g., it is 4.77 for the 
"Pure_C_Xtal" case. 

correction given explicitly in Table 4 of iRicher et al.l (120081) 
and shown in the top panel of our Figure 3. Finally, we nor- 
malize the resulting curves by minimizing their root-mean- 
square residuals in the neighborhood of the peak, between an 
F814W of 25.1 and 27.5. 

In Figure 4 we show the LF of our fiducial sequence 
("Pure_C_Xtal"), that of a sequence with crystallization ar- 
tificially suppressed ("Pure_C_No_Xtal"), and that of one 
including crystallization but artificially excluding the latent 
heat of crystallization ("Pure_C_Xtal_No_LH"); all are plot- 
ted over a histogram of the observed LF. The no crystallization 
and no latent heat sequences show evidence of a bump due to 
convective coupling around F% 14W ~ 26, but do not continue 
to rise through the observed maximum; this is clearly incon- 
sistent with the data, and no adjustment of mass or internal 
composition can bring them into good agreement. In terms 
of x 2 , for average observational errors of ~ 5.5 stars/bin in 
the neighborhood of the peak, we have \ 2 = 0.75 for the crys- 
tallizing sequence and \ 2 ^2.2 for the no crystallization and 
no latent heat sequences, a nearly three-fold increase in x 2 - 
Thus, the sequence with crystallization provides a much bet- 
ter match to the data. 

3.3. Constraining crystallization, phase separation, and core 
composition 

In the years since Ivan Hornl (1 19681) the realization that 
the cores of normal mass WD stars should consist of 
a mixture of carbon and oxygen implied that crystal- 
lization may also release energy result ing from phase 
separation of the carbon and oxygen dStevensonl 119771; 
i Barrat. Hansen. & Mpchkoyitchlll988t [S egretain & Chabrier 
ll993HSegretain et aljl994tllsern et alj2000h . This occurs be- 
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FIG. 5. — The same as Figure [4] but for DA model sequences with 
uniform, 50:50 carbon/oxygen cores: crystallization only ("CO_Xtal"), 
crystallization and phase separation ("CO_Xtal_PS"), and no crystalliza- 
tion ("CO_No_Xtal"). In addition, we show a pure oxygen DA sequence 
("pure_0_Xtal") and a pure carbon DB sequence ("pure_C_DB_Xtal"). All 
models have 0.5Mq . 

cause when a carbon/oxygen mixture crystallizes the oxygen 
content of the solid should be enhanced. Since WDs crystal- 
lize from the center outwards, this leads to a net transport of 
oxygen inward and carbon outward, and because oxygen is 
slightly heavier than carbon this differentiation releases grav- 
itational energy. 

We hav e included this en e rgy in our models as de- 
scribed in Montgo mery et al.l d 1999b . For these compu- 
tations we have assumed the carbon and oxygen abun- 
dances are equal throughout the core. This underesti- 
mates the oxygen abundance compared to th at predicted 
by sta ndard stellar evolution calculations (e.g., ISalaris et alJ 
119971) . bu t the remaining uncertainty in the C(q,7)Q reac- 
tion r ate dMetcalfe. Salaris. & Wingej|2002l: lAssuncao et alJ 
2006) lead to a degree of uncertainty in the C/O ratio and 
profile. In Figure 5 we show several LFs: crystalliza- 
tion only ("CO_Xtal"), crystallization with phase separation 
("CO_Xtal_PS"), no crystallization ("CO_No_Xtal"), pure 
oxygen core with crystallization ("Pure_0_Xtal"), and a pure 
carbon DB sequence ("Pure_C_DB_Xtal"). As is readily ap- 
parent, all of these sequences have a peak which is too bright 
by at least 0.5 mag in F 8 14W. 

These results have several interesting possible interpreta- 
tions. First, the results seem to suggest that the oxygen con- 
tent of these stars is relatively small or zero, since it is the 
higher crystallization temperature of oxygen which shifts the 
peak in the LF to smaller magnitudes. The only way to 
accommodate more oxygen would be to hav e lower mass 
mode ls, in conflict with the distance modulus (lHansen et alJ 
120071) . The colors (e.g., Figure 2) also make it difficult to ap- 
peal to lower masses with higher oxygen abundances. Ad- 
ditionally, for plausible IFMRs the main sequence lifetime 
for single stars becomes more problematic with lower WD 
masses even with the IFM R dependency on metallicity of 
iMeng. Chen. & Hani d2008l) . Thus the constraint on the inte- 



rior oxygen abundance becomes stronger. Taken at face value, 
these results indicate that the carbon-to-oxygen ratio is much 
greater than 1, and we will be able to make a more quantitative 
statement from our future more complete Bayesian statistical 
analysis (in preparation). 

Second, as shown in Figure 4, the data are consistent and 
well-fit by carbon core models with crystallization and the re- 
lease of latent h eat, but not by mo dels without. This confirms 
the prediction of Ivan Hornl (Il968h that crystallization is a first 
order phase transition and releases the latent heat of crystal- 
lization. Were it not so, crystallization would leave no sharp 
peak at this magnitude in the observed LF. This impacts our 
understanding of solid-state physics at extremely high den- 
sity: it is the first empirical confirmation of the release of 
latent heat during crystallization - an important theory that 
has a large impa ct on WD ages, as has been pointed out by 
Ivan Horn! (1 19681) and many authors in the intervening years. 

Third, it is possible a priori that a significant fraction of the 
observed WDs may be DBs. The mismatch of observed LF 
and the DB sequence in Figure 5 essentially eliminates pure 
He atmosp heres as a significan t component of the sample (as 
shown by lHansen et al.l 120071) . but not models that become 
mixed (H/He) as they cool; we explore this possibility in a 
forthcoming paper. 

Fourth, the "best-fit" fiducial sequence in Figure 4 begins 
crystallizing near the value of r xta i = £couiomb/£'therm a i ~ 170. 
If the actual value is higher (lower) then crystallization will 
occur at lower (higher) luminosities. Higher values allow fits 
with larger amounts of o xygen in their cores. In addition, 
Potekhin & Chabrier (2000) show that even for a pure compo- 
sition theoretical uncertainties in the polarization of the elec- 
tron Fermi gas and quantum effects in the liquid and solid 
phase can alter the value of r xta i- In future analyses, an accu- 
rate determination of the mass, distance, and reddening will 
lead to an accurate determination of r xta i and the core compo- 
sition. 

Finally, we note that the central density and temperature as- 
sociated with a particular value of F814W through the model 
atmospheres is sensitive only to the mass-radius relationship 
set by the degenerate electron pressure support. This is very 
insensitive to the C/O relative abundances — these produce 
differences of SFRIAW < 0.05. This implies that the value of 
r in the center, at the peak of the LF for example, is sensi- 
tive only to the interior composition. Therefore we conclude 
that the onset of crystallization is determined by the particular 
mixture and the value of V for that mixture. Comparison of 
the theoretical models and the data promise to provide impor- 
tant measures of the onset and development of crystallization. 

4. DISCUSSION, SUMMARY, AND FUTURES OF EXPLORING 
WHITE DWARF PHYSICS WITH CMDS 

Although we are not focused on uncertainties, it is reason- 
able to examine how changing the distance modulus might 
affect the results. Put another way, how much does the dis- 
tance have to change to reproduce the peak of the observed 
LF with oxygen crystallization rather than carbon? The an- 
swer is contained in Figure 5. Here we see that to make the 
location of the peak of the LF consistent with oxygen crys- 
tallization we have to lower the distance modulus by a little 
more than 0.5 magnitude s — this possibility is excluded by 
the main sequence fitting (R icher et alJl2008T) . 

We have shown that simultaneously fitting the main se- 
quence and the WDs in a cluster gives the best possible 
constraint on distance, metallicity and reddening corrections. 
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Physically realistic atmosphere calculations then allow us to 
place evolutionary tracks in the CMD. The number distribu- 
tion of stars contains important information on the internal 
physics of the WD stars. This allows us to explore the physics 
of crystallization. We present evidence that the data is most 
consistent with a first order phase transition , releasing laten t 
heat during crystallization, as proposed bv lvan Hornl 1l968). 
The current data places constraints on the onset of crystalliza- 
tion, the central carbon/oxygen abundance, and the compo- 
sition of the envelope at the degeneracy boundary. We will 
improve these constraints with a more complete Bayesian sta- 
tistical analysis in the near future. This work also points to the 
importance of forthcoming data on additional clusters as well 
as increasing the sample of stars through more HST fields on 
this cluster. This work also underscores the essential nature 
of more proper motion data to get the most information out of 
these kinds of studies. 

Pulsations may also allow an asteroseismological determi- 
nation of th e crystallized mass fraction for massive pu lsators, 
as shown by Metcalfe. Montgomery, & Kanaan (2004) for the 
DAV BPM 37093, although this claim has been challenged by 



iBrassard & Fontaine! d2005l) . While certainly important, we 
note that asteroseismological analyses do not probe the latent 
heat of crystallization; this quantity is accessible only through 
the WDLF, as demonstrated in this paper. To this end, we 
eagerly anticipate the forthcoming HST observations of this 
cluster. These will provide a proper-motion screened sample 
over a larger area of the cluster and to fainter magnitudes, pro- 
viding an exacting test of the ideas put forth in this paper. 
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